PARAFERMION STATISTICS AND THE APPLICATION TO NON-ABELIAN 

QUANTUM HALL STATES 

EDDY ARDONNE 

Abstract. The (exclusion) statistics of parafermions is used to study degeneracies of quasiholes 
over the paired (or in general clustered) quantum Hall states. Focus is on the TL^ and su(3)fc / u(l)'^ 
parafermions, which are used in the description of spin-polarized and spin-singled clustered 
quantum Hall states. 



1. Introduction 

In the last decade, low-dimensional systems in which the fundamental (quasi) particles do not 
satisfy bose or fermi statistics have received a lot of attention. Among the most famous examples 
are the fractional quantum Hall (fqH) systems. The quantum Hall states at simple filling fractions 
V = where M is an odd integer, are understood in terms of the famous Laughlin states 
0. Quasiparticles (or quasiholes) over these states carry fractional charge and satisfy fractional 
statistics. 

In the last few years, generalizations of these states have been under investigation. Among these 
generalizations are so-called paired (or more general, clustered) quantum Hall states. The notion 
of clustering will be explained in the next section. In the definition of these states, parafermions 
play a predominant role. Ultimately, it is the presence of these parafermions which cause the 
quasiholes to have peculiar statistics properties, namely, they obey 'non-abelian' statistics (see, for 
an early reference, Q). 

The simplest of the paired quantum Hall states is the pfaffian state proposed by G. Moore and 
N. Read |^ ; it is believed that this states describes the quantum Hall effect at filling fraction v 
This quantum Hall effect is special in the sense that it is the only quantum Hall effect observed 
at a filling fraction with an even denominator (in single layer samples) Q . For a recent review 
on this subject, see |5]. Clustered analogues of the paired pfaffian state were proposed by N. Read 
and E. Rezayi (RR) 1 1 . The states mentioned above are all spin-polarized; spin-singlet analogues 
were proposed in Q . Recently, the bosonic versions of the RR states were shown to be relevant in 
the context of rotating Bose-Einstein condensates js] . 

The statistics properties of the parafermion fields will be investigated in this paper, with the 
intention of obtaining closed form expressions for the ground state degeneracy of clustered quantum 
Hall states in the presence of quasihole excitations, as described in This state counting 

problem is interesting for the following reasons. The clustered quantum Hall states can be seen 
as ground states of a hamiltonian with an (ultra local) interaction between the electrons. Finding 
the ground state degeneracy of this hamiltonian can be done in a conformal field theory (CFT) 
approach, relying heavily on the statistics properties of the parafermionic fields. Another approach 
is by numerically diagonalizing the interaction hamiltonian for a small number of electrons. This 
method can serve as a check on the analytical results of the first approach. Thus, the quasihole 
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degeneracies of a system of interacting electrons can be understood in terms of parafermionic 
statistics! 

In the context of the the spin-polarized states of Read and Rezayi, the (or su(2)fe/u(l)) 
parafermions are the relevant parafermions. For the 'non-abelian spin-singlet' (NASS) states of 

0, the relevant parafermions are the parafermions related to su(3)fc/u(l)^ (see for a discussion 
on general parafermion CFTs). 

The plan of the paper is as follows. We first recall in which way (clustered) quantum Hall 
states can be constructed in CFTs (section ||) . In section || we will shortly indicate the setup of 
numerical diagonalization studies, because we need to adapt the calculations to the setup in which 
these studies are done. The general structure of the counting formulas will be indicated in section 
^. It will become clear that the degeneracy consists of an intrinsic and an orbital part, which 
need to be combined in the right way. The intrinsic degeneracy factors need to be split to make 
this possible. The remainder of the paper is devoted to this task, as no explicit expressions for 
these 'split degeneracies' were known for general level k > 2 for the states under consideration 
in this paper. We will explain the procedure to obtain these expressions using the su(3)2/u(l)^ 
parafermions of the NASS state at level A: = 2 as an example (see fl^). The first step is to find 
a basis for the (chiral) spectrum of the parafermion CFT. Using this basis, recursion relations for 
truncated characters will be derived (section ^). These recursion relations can be solved using the 
results of section^, providing expressions for the truncated characters. From the explicit truncated 
characters, the 'split degeneracies' can be extracted. Finally the counting formula for the paired 
spin-singlet states is obtained in section 0, filling in some of the details of the discussion in lll^ . In 
section ^, counting formulas for the RR states at general level k are obtained, while section gdeals 
with the counting formulas for the general k NASS states. Section |l^ is devoted to discussions. 
In this paper, we do not go into the details of the numerical analysis, but refer to the papers 
p2| , H, H, For all the cases checked, the results of the numerical diagonalization studies are 
exactly reproduced by the counting formulas. 

2. Construction of clustered quantum Hall states 

Clustered quantum Hall states are constructed as correlators in certain CFTs. But before we 
come to this construction, let us first explain the nomenclature 'clustered states'. Clustered states 
(with order k clustering) have the following form 

^¥iU^})^<osiU^})Y[('^~^,f' , (1) 

where $bos ^ bosonic factor (symmetric under the exchange of any two coordinates) with the 
property 

I $L(^i = ^2 = --- = ^z)^0 i<k 

\ ^t,{zi =Z2 = --- = Z^)=0 i>k , ^ ' 

1. e. the wave function vanishes as soon as fc -I- 1 or more particles are at the same location. Note 
that we omitted the exponential factors of the wave functions. The factor ni<j('^* ~ ^j)^^ is the 
familiar Laughlin factor. The fermionic quantum Hall states all have M odd, while the bosonic 
versions would have M even. Though changing M changes many properties of the states, such as 
the filling fraction, the quasihole charge and statistics, the quasihole degeneracies are unaffected, 
and the counting formulas presented in this paper hold for any M. Note that this degeneracy is 
not to be confused with the 'torus degeneracy', which in fact does depend on M. 

Following the pioneering paper quantum Hall states can, under certain restrictions, be 
defined as correlators in CFTs. For the clustered quantum Hall states, CFTs with affine Lie algebra 
symmetry are used. The operators creating the electrons in general consist of a parafermion field 
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and a vertex operator of a set of r chiral boson fields (r is the rank of the afhne Lie algebra). These 
boson fields take care of the charge, spin and possibly other quantum numbers associated to the 
particles. The electron creation operators take the form 

V,,{z,) = M^.i):e'f-'^{z,):, (3) 

where ipa is a parafermion field and <p is a set of free boson fields. The quantum Hall states can 
be now be defined as follows 

hm z^(youK=i.2---K=i,iV :e-''b'^(zo,) :). (4) 

Zoo— 

The parameter b in the background charge needs to be chosen in such a way that the overall 
correlator is charge neutral, a is chosen such that the effect of the background charge does not go 
to zero in the limit Zqo oo, whilst keeping the result finite. Moreover, the (unique) fusion of the 
parafermion fields ip must result in the identity operator 1 in order to get a non-zero correlator. 
This in general puts a restriction on the number of electrons N. In this paper, we will not give 
explicit results for the correlators (and thereby the wave functions '^^)- These wave functions 
and other properties of the ground state, such as the filling fraction, can be found in for the 
Read-Rezayi states and in for the spin-singlet analogues proposed in . A different form of the 
Read-Rezayi wave functions can be found in [T^ . Also, wave functions for states with quasiholes 
(see below) are presented there. 

States with excitations (quasiholes) present are also defined in terms of correlators in the same 
CFT as the parent state, by inserting the corresponding quasihole operator in the correlator. Cer- 
tain properties of these quasiholes can be studied via the corresponding correlators. A constraint 
on these operators follows from the condition that the wave function has to be analytical in the 
electron coordinates (the lowest Landau level condition). This implies that the quasihole operators 
need to be local with respect to the electron operators. This constrains the quasihole operators to 
be of the following form 

V^p{w) = : e'l^'-'^iw) : , (5) 

where is a 'spin field' of the parafermion CFT. Note that the nomenclature spin field does 
not refer to the electron spin. Again, one has to insert a background charge to enforce charge 
neutrality. Also, upon fusing the parafermion and spin fields, one has to fuse to the identity 
operator 1 in the last step, to obtain a non-zero correlator. This constrains the possible particles 
in the correlator. The fusion of the spin fields is not unique: in general there is more than one 
fusion channel. This implies that a correlator with several quasihole operators in general stands 
for more than one quantum Hall state. In other words, the clustered quantum Hall states with 
quasiholes are degenerate. It is this degeneracy which lies at the heart of the non-abelian statistics. 

One immediate question one can ask is how many states does one describe with such correlators? 
In fact, this question can be answered in two independent ways. The first one is via a numerical 
diagonalization study of interacting electrons on the sphere, in the presence of a magnetic field. 
In this paper, we will follow the second approach, which is analytical, and uses the conformal field 
theory of the underlying parafermions. But before we come to this point, we have to spend some 
words on the numerical approach as well, in order to be able to adapt the analytical approach to 
the numerical setup. This is the subject of the next section; for a more details, we refer to ||ic[| . 

3. The setup of the numerical studies 

Though we will not describe numerical diagonalization studies in depth in this paper, it is nec- 
essary to point out briefiy in which setup they are done, because we need to adapt our calculations 
to be able to compare results. The numerical diagonalization is most easily done on the sphere. 
The interaction between the electrons is chosen such that the clustered state under investigation is 
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the unique ground state (in the absence of quasihole excitations). Note that this interaction is an 
ultra local, many-body interaction, rather different from the long range Coulomb interaction. To 
'tune' to the right filling fraction, a specific number of flux quanta need to penetrate the sphere. 
States with quasiholes can be studied by increasing the number of flux quanta (but keeping all 
the other parameters the same); this results in the creation of quasiholes, as can be seen from 
the Laughlin gauge argument. The number of flux quanta needed for a state on the sphere with 
quasiholes is given by 

N^ = ^N-S + AN^ , (6) 

where N is the total number of electrons, and AA'^ the number of excess flux quanta, needed for 
the creation of the quasiholes. S is an integer constant depending on the state under investigation. 
Also, the number of quasiholes which are created by increasing the number of flux quanta by one 
depends on the state under investigation. For the spin-polarized RR states, this relation is given 
by n = kANcj,, where n is the number of quasiholes. For the spin-singlet analogues, we have 
n = n-^ + ni = 2kAN^. 

For the clustered quantum Hall states with quasiholes present, the ground state is degenerate 
(for the ultra local interaction). The degeneracy consists of two parts. First of all, there is an 
orbital degeneracy, which is caused by the fact that in this setup, the quasiholes are non-local. 
This orbital degeneracy is not speciflc for clustered states; it is also present for the (unpaired) 
Laughlin states. For a system in which the quasiholes are localized, this degeneracy would not be 
present. Secondly, there is an intrinsic degeneracy, which stems from the non-trivial fusion rules 
of the spin fields, needed to create quasihole excitations. This source of degeneracy is special for 
the clustered states. In this paper, we will focus on this intrinsic degeneracy and obtain analytical 
expressions, which allow the combination with the orbital degeneracy factors. This provides us 
with explicit expressions for the degeneracy of the ground states, in the presence of quasiholes. 

As spin and angular momentum are good quantum numbers, all the states obtained from the 
numerical diagonalization fall into spin and angular momentum multiplets. The structure of the 
counting formulas is such that also the multiplet structure can be extracted. 

4. Degeneracy factors and counting formulas 

The intrinsic degeneracy is caused by the non-trivial fusion rules of the spin fields. As an 
example, we will use the spin fields of the 5u(3)2/u(l)^ parafermionic CFT. The fields and their 
fusion rules in this theory can be determined according to the methods of pT[ | and are summarized 
in table |l|. We use the notation introduced in The parafermion fields are denoted by tp, and 
all have conformal dimension = i. In particular, ipi,tfj2,4'i2 correspond to the roots ai,—a2 
and ai+a2 of sm(3), respectively. The spin fields (T|,(T|,(T3 and p are related to the weights of sm(3) 
and their conformal dimensions are given by A^. = and Ap = |. The fusion of an arbitrary 
number of a-\^i fields can be depicted in a Bratteli diagram (see also [^). Each arrow in the 




Figure 1. Bratteh diagram for the spin fields of su(3)2/u(l)^. 
diagram in figure |l| stands for either a CT| or ai field. The arrow starts at a certain field which 
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Table 1. Fusion rules of the parafemiion and spin fields associated to the 
parafermion theory su(3)2/u(l)^ introduced by Gepner [ pi] . 



can only be one of the fields on the left of the diagram at the same height. This last field is fused 
with the one corresponding to the arrow, while the arrow points at a field present in this fusion. 
As an example, the arrows starting at the * are encoding the fusion rules p x (J^{i) = V'2(i) + f T(i) 
and cr3 X cr|(|) = 01(2) + o'i(T)- One checks that the diagram is in accordance with the first two 
columns of Table ^ 

From figure |l|, one immediately reads off that in general there is more than one fusion path of 
spin fields with leads to the identity (possibly the identity is reached only after the fusion with 
the parafermion fields 0i_2 of the electron operators). It is easily seen that the number of fusion 
channels starting from and terminating at 1 while rt| and tT| spin fields are fused is given by 

dn-r,n; = ^(?T-T +ni-2), (7) 

where J-{n) is the n-th Fibonacci number, defined by J-{n) = J-{n — 1) + J-{n — 2) with the 
initial conditions ^{0) = 1 and ^(1) = 1. Next to this intrinsic degeneracy, there is an orbital 
degeneracy. These orbital degeneracy factors can be found in ^ for the states discussed in this 
paper. These factors have the general form 

The product is over the types of quasiholes, while the numbers Fi are interpreted as the number 
of 'unclustered' particles in the state. In the correlators, these correspond to the fundamental 
parafermions ipi. For each fusion path, these numbers can be different, implying that we have 
to split the intrinsic degeneracy according to these numbers. We denote these 'split degeneracy 
factors' by Explicitly, we have { ]^ }a; and { ^[ }k for the RR and NASS states, respectively. 

Using the notation above, the counting formula for the clustered spin-singlet quantum Hall 
states take the following form 

where the prime on the summation indicates the presence of constraints (see below eq. (^l|)). The 
equivalent counting formula for the Read-Rezayi states is given in eq. (^3|). 

The new result of this paper are the explicit split degeneracy factors for the and 5u(3)fc/u(l)'^ 
parafermions at level k > 2 (for fc = 2, these factors can be found in [|2| and Q respectively). 
Previously, these factors for the Zfc parafermions (fc > 2) could only be characterized via recursion 
relations, see Note that the results in this paper are easily extended to the more general 

su(iV)fc/u(l)^~^ parafermions. 
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We will now briefly outline in which way the split degeneracy factors are obtained. The starting 
point is the character of the parafermionic CFT. The symbols {}k can be extracted from finitized 
forms of these characters |^ (see also ) . Recursion relations for these finitized characters can be 
derived from an explicit basis of the parafermionic CFT. These recursion relations will be written 
in a way that allows for an explicit solution, from which the symbols {}k can be extracted. In the 
sections we will demonstrate this for the level k = 2 spin-singlet states of . 

5. A BASIS FOR THE 5u(3)2/u(l)^ PARAFERMION THEORY 

In this section, we briefly describe how an exphcit basis for the chiral spectrum of the su(3)2/u(l)^ 
parafermion CFT is formed. The starting point is the chiral character for the parafermions in the 
su(3)2/u(l)^ conformal field theory. This character can be written in the form of a 'Universal 
Chiral Partition Function' (UCPF) sec, for instance, ||l6|, |l^. This character reads 

ch(a;i,a;2;g,fc = 2) = 2^ — — — xl^x^'' ■ (10) 

In this character, Xi = e'^'^' are fugacities of the particles, and q = e^^ {(3 is the inverse tempera- 
ture). {q)a is defined by {q)a = nfc=i(l ~ for a > and {q)o — 1. 
The bilinear form in the exponent of q is described by the matrix 



A 7 

2 ^ 



(11) 



The same matrix also describes the exclusion statistics of these parafermions. More information 
on the relation between exclusion statistics and the UCPF can be found in . 

A basis for a CFT can be thought of as a set of states spanning the chiral Hilbert space. This 
set of states can be written as a (set of) vacuum state(s), on which creation operators act. The 
parafermions V'i,2(^) in the su(3)2/u(l)^ theory can be expanded in modes as 



^1,2 (^) = E • (12) 



As usual, the modes tpm with negative index are the creation operators while the modes with 
positive index annihilate the vacuum 

V'm|0)=0 m>0. (13) 

The set of states 

V-".„Cr„^,---^"^.jo) (14) 

is overcomplete, because of the (generalized) commutation rules of the parafermions. In the fol- 
lowing, we will point out which restrictions on the indices Si will remove the 'overcompleteness'. 
In doing so, we will follow the exclusion interpretation of the K-matrix as closely as possible and 
concentrate on the lowest possible 'energy' (given by La = J^i ^i) ^ certain number of applied 
fields first. The ordering of the modes '0^'^ is such that we apply the -0^ modes first. From ( p^ ) it 
follows that the simplest non-trivial state is 

^ii/2|0) . (15) 

Interpreting the matrix ( pT| ) as the exclusion statistics matrix, the minimal spacing between two 
il>^ modes is 1, thus the state with two ^p^'s acting on the vacuum with minimal energy is 

^-3/2V'il/2|0) . (16) 
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The extension to rii ip^ modes is simple 

^-(2,n-l)/2---i^-V2^\/2\0) ■ (17) 

Note that if this were the whole story, we would describe the (free) Majorana fermion. The spacing 
between ip'^ modes is the same as for the ip'^ modes. However, if one acts with ip'^ on a state in 
which modes are already present, one has to take into account the mutual statistics between 
and -0^ modes, which is —1/2, according to dTl]). Thus the energies of the ip'^ modes have an 
extra shift of — ni/2, resulting in the following states (with minimal energy) 

'>P-(2n2-l-,n)/2 ■ ■ ■ V'-(3-„i)/2'/'-(l-ni)/2^i(2rn-l)/2 ' ' ' ^'^3/2^-^1/210) . (18) 

The (dimensionless) energy associated to this state is ^±+'^2-^1^2 ^ ppgf-jggiy exponent of q in 
the character (p^. To obtain all the possible states, one has to allow states with higher energies 
as well. As usual jlj], the energies of all the modes can have integer shifts, under the restriction 
that modes acting on a state have larger energies than the modes of the same type which have 
been applied earlier. This results in the following set of states 

■^-(2n2-l-ni)/2-t„2 ' ' ' V'^(3_„i)/2-t2 ^-(1-ni )/2-ti V^-(2ni -l)/2-s„^ ' ' ' "^^3/2-32 ^-1/2-si 1 0) ' 

(19) 

with s„^ > . . . > S2 > si > and tn^ > ■ ■ ■ > t2 > h > {si,tj G N). 

Up to now, wc used the special ordering of applying modes to the vacuum, namely, all the 'ip^ 
modes first. This is in fact enough to span the whole chiral spectrum, as can be seen if we perform 
the trace over all basis states. More or less by construction, we obtain the character (|l^). However, 
we also can allow a general ordering of the modes. As an example, we take the following state 

^-o^-i/2\Q) ■ (20) 

The energy of the tp^ mode is zero because it gets an extra shift of —1/2 due to the presence of 
the ip^ mode. In spanning the whole chiral spectrum, we can also choose to use the state, with the 
order of the modes changed 

^-o^-i/2|0) . (21) 

In this case, the ip'^ mode gets an extra shift of —1/2, because of the presence of the ^'^ mode. 
Thus, the Lq value is the same for both states. In general, changing the order of two neighbouring 
ip^ and ip'^ modes does not change the Lq value if the extra shifts are changed in the appropriate 
way. The extra shift of a field is given by —1/2 times the number of preceding modes of the other 
type. In general, two states related by a reordering of modes are different, but we can use either of 
them (but not both) to span the chiral spectrum. Note that the rules of the spacing between the 
various fields is in accordance with the (exclusion) statistics interpretation of the matrix K. The 
character (|l^) is obtained by taking the trace over all the states in the basis ( |l9| ) 

ch(xi , X2 ; g) = Tr x"^ x"^^ . (22) 

We can now define the finitized characters needed in the derivation of the symbols {\k by using 
the basis described above. These finitized characters are polynomials which will be denoted by 
Y{i^m)- These polynomials are traces over the basis (|9|), but restricted to the states in which the 
energy of the modes of the ipi (-02) fields are smaller or equal to I (m). Though the total energy 
of a state does not depend on the ordering of the modes, the energies of the individual modes do 
depend on the ordering, as can be seen by comparing the states ( pO| ) and (|2^). By restricting the 
trace over states in which the labels of the modes are bounded, we must include a state if there is 
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at least one ordering in which aU the modes satisfy the bounds imposed. Note that there may be 
other orderings, in which these bounds are not satisfied. We write the finitized characters as 

Y(i,^){xi,X2;q)^ ^-Tr' x^^"'?''" • (23) 

<l,<m 

The prime on the trace denotes an important restriction on the number of modes (denoted by ni 
and ^2) present in the states. These numbers must satisfy ni — 21 (mod 2) and 712 — 2m (mod 2). 
This restriction takes into account that after fusing the spin fields, one ends up in the right sector, 
which can be l,7/'i,7/'2 or depending on the number of spin-up and down electrons. This is 
necessary, because after fusing the spin fields and the parafermion fields of the electron operators, 
on has to end with the identity 1, to obtain a non-zero correlator. 

The finitized characters ^(/.m) can be written in terms of recursion relations of the following 
form 

y{i,m) = Y^l^„,_^■^+X2q"'~^Y(l^l„^_^). (24) 

Note that the recursion relations above are stated in terms of the energy labels of the modes. The 
aim we have is finding the number of possible states when a certain number of extra flux is added. 
We therefore need to make a change to labels which depend on the additional flux. In fact, we will 
use the number of particles (given by n-\ and in this case) created by this flux as labels for the 
finitized partition functions. Explicitly, we have I = ^ and m = In terms of the number of 
created quasiholes, the recursion relations become (compare [p^) 

^(n-r,n;) = ^(n-r,n^-2) + 2;2g ^ Y(„^ + i^„^_2) . (25) 

The initial conditions for these recursion relations look as follows 

Y{i,Q) = Y{o,i) = , 

^(0,0) = ^(2,0) = ^(0,2) = 1 J 

Y(x.i) = q^xiX2 ■ (26) 

The finitized characters are completely described by ( psf ) and (^6|). In the next section, we will 
solve these recursion relations and thereby provide explicit expressions for the finitized characters. 



6. Recursion relations and solutions 



The recursion relations of the previous section can be solved explicitly; we will follow the 
approach of |Q. The key observation is that the recursion relations can be matched to general 
recursion relations, which are solved in terms of finitizations of universal chiral partition functions. 
Consider the following polynomials P 



P^iz;q)=J2(l[<")l 



^ m ■ K ■ m + Q ■ m 



n 



(L + (I - K) • m + u) 



(27) 



In this equation, I is the identity matrix, K the statistics matrix and [ ^ ] the g-deformed binomial 
(q-binomial) 




(«)a 



a,b eN; b < a 
otherwise 



(28) 
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Note that we defined the g-binomial to be non-zero only if both entries are integers greater or 
equal to zero, to avoid additional constraints on the sums in the counting formulas. 
From the definition of the g-binomials, the following identity is easily derived 

Replacing the i'th q-binomial factor in ( ^7| ) by the right hand side of (|29|), one finds the following 
recursion relations 



a 




'a - I 




a — 
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b 




b 




b- 


1 



PL(z;q) = PL-e,(z;q) + z.g-^^''+'^'+"'+^'PL-K.e,(z;'z) • (30) 

The vector represents a unit vector in the z'th direction. We will use the equivalence between 
( p7| ) and ( |30| ) frequently, because the recursion relations we encounter in deriving the counting 
formulas are all of type (pO|). Of course, upon deriving polynomials from recursion relations, one 
has to take the initial conditions into account. For the counting we need to know the finitizations 
of the character formulas, and these can be written in the form (|2^). Thus, when we solve recursion 
relations by polynomials of the form (|2^), the proper initial conditions are automatically taken 
into account. 

We start by applying the above to the recursion relations (|2^), resulting in the following ex- 
pressions for the truncated characters Y(„^^„^) 

-ni -t-a- 



y{nT,rn)ixi,X2;q) 



{a^+b^-ab)/2 a b 



a,b 



rif+b - 
2 



(31) 



This result will be needed for the final counting formula, which we give in the next section. 

7. A COUNTING FORMULA FOR THE NASS STATE AT fc = 2 

From the truncated characters of the previous section, we can obtain the symbols {}2, needed 
in the counting formula eq. (||). In fact, the symbols {}2 are obtained by taking the limit (7 — > 1 

in eq. (|3l|) (see, for instance, |0, |lOt) 



of the coefficient of xf^ x^" 



F2 



In this limit, the q-binomials in 



l|) become 'ordinary' binomials and we find 

T+-F2 \ /rii+£i~ 



F2 



2 

Fi 



2 

F2 



(32) 



(33) 



The fact that the finitized characters indeed provide the symbols {} is rather non-trivial. This 
connection was first proposed in Some (restricted) 'solid on solid' (SOS) models (see, for in- 
stance 1^) can be mapped to the Bratteli diagrams of the spin fields of the quasiholes. Recursion 
relations for the partition functions for these models (at finite size) are in general equivalent to 
recursion relations for finitized characters in certain CFTs. In the case at hand, the correspond- 
ing CFT is the parafermion CFT. This provides a link between the Bratteli diagrams and the 
parafermion theories. As a check, on can calculate the number of fusion paths for the spin fields by 
summing over the symbols {} and compare to the result obtained from the diagram itself. In this 
specific case, the equivalence follows from the structure of the recursion relations (see for instance 
pO|), giving rise to the identity 

-Ti ,-r2 
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The prime on the summation denotes the constraints -Fi = (mod 2) and F2 = (mod 2). At 
the level of the wave functions, the degeneracy is due to the presence of particles which do not 
belong to a cluster any more. At the level of correlators, these unclustered particles correspond to 
parafermions V'^ and ijP' , which act as 'cluster breakers'. In the case of the pfaffian state, this was 
made explicit in Q . 

The counting formula for the NASS state at fc = 2 is obtained by inserting the symbol (^) in 
the general counting formula (^ 

-'''T,4-i"T, 1.-^^1.2 

where the prime on the sum indicates the constraints + Ni = N, n| + = 4AiV0 and 

We will now comment on the spin and angular momentum multiplct structure. As an example, 
we will write out the polynomial ^(7,1) 

^(7,1) = il^ +qi + q^)xiX2 + [q^ + 2qi + 2q^ + 2q^ + q^)xlx2 + q^ xlxl . (36) 

After multiplying the coefficient of xf^x^^ with (in general) g-("T^i+"-4--P'2)/4^ obtains a sum 
of terms of the form g'*" , which together form a collection of angular momentum multiplets with 
quantum numbers Iz- For instance, the coefhcient of x\x2 in eq. ( p6[ ) encodes two multiplets, 
namely L = 2 and L = 1. 

An alternative way to obtain these results is by associating angular momentum multiplets to 
the binomials in eq. (|3g). The binomials forming the symbols {}2 need to be interpreted as the 

number of ways on can put / fermions in a boxes, which are labeled with Iz = — 1 — + 
1, . . . , -^^Y^ angular momentum quantum numbers. Each way of putting the / fermions in a boxes 
has an associated with it. Together, these quantum numbers fall into angular momentum 
multiplets. In this way, angular momentum multiplets can be associated with the binomials. The 
angular momentum multiplets of the various binomials in the counting formula need to be added 
in the usual way. 

Though the parafermion theory does not have a proper SU (2) spin symmetry, one can associate 
spin quantum numbers to every state by taking Sz = — -■ Combining the spin and angular 
momentum, one finds that all the states fall into spin and angular momentum multiplets. 

The numerical diagonalization studies for the NASS state at level k — 2 are described in [l0[| . 
It is very gratifying to see that the counting formula eq. (|3^) does in fact exactly reproduce the 
quasihole degeneracies, as well as the multiplet structure. 

In order to find the counting results for the spin-singlet states at general level-A;, we first take a 
closer look at the counting of the Read-Rezayi states, which was in fact done in Those results 
however, were stated in terms of recursion relations which are difficult to solve. The advantage of 
the recursion relations presented in the next section is that they can easily be solved in terms of 
(q-deformed) binomials, and thus provide explicit expressions for the symbols {}k- 



8. Counting formulas for the Read-Rezayi states 

The derivation of the counting formulas for the RR states goes along the same lines as the 
derivation for the NASS k=2 states as explained in the previous sections. Therefore, we do not go 
into full detail, but concentrate on the parts which need more explanation. 
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We start with the character of the su(2)fc/u(l) parafermionic theory (see |g2[ 
obtained from ||, H 



ch(a;; q, k) 



which can be 



(37) 



where a — (ai, . . . , ak-i) and C^^i = 2Aj.j';-|^, 
nents, these matrices are given by 

(Afc_i)jj — 2(5jj — <5|i_j|_i 
(K-ikj = min(i,j) - ^ 



_i being the Cartan matrix of su(k). In compo- 



i, j 1, . . . ,/j - 1 



(38) 



In fact, Cfe_i is the K-matrix for the parafermions, see |17|] . The parafermions in this theory 
are t/iqj V^ij • ■ • V'fc-i (V'o is the identity 1 and the labels are defined modulo k). By writing x*"' in 
the character (^^, we take care of the fact that the fugacity of species i is i times the fugacity 
of the first type of particle. In fact, the i'th species can be thought of as a 'composite' of i 
particles of species 1. This point of view is supported by the fusion rules for these parafermions 
ipi X = V'p+ij with p=l,...,fc— 1. This structure is also present in the K-matrix structure 
describing the Read-Rezayi states |2^, |l^ . 

A basis for the chiral spectrum can be constructed in the same way as described in section ^. 
The shifts in modes between the various fields are given by the elements of the matrix K. We will 
now proceed by directly giving the corresponding recursion relations 



Yi{x;q,k)=Yi^ 



X q 



Ji-- 



"11- 



(39) 



The factor 



-J, — is the conformal dimension of the i'th parafermion in the Z^-parafermion theory. 
These recursion relations are solved by the following polynomials 



Y,(x;g,fc)-^g^(a-C 



k-l 

n 



(l+(Ife_i-Cfc_i)-a), 



(40) 



where Ifc_i denotes the (k — l)-dimensional unit matrix. To obtain the counting results, we have 
to specify the truncation parameters h. As in the NASS case for fc = 2, we will do this in terms 
of the number of particles created by the extra flux, given by n = kAN^ for the states under 
consideration. Because the chemical potential of species j is i times the chemical potential of 
species 1, the truncation parameter li has to be scaled with a factor i with respect to h (see, for 
instance |19 ), which is found to be Zi = ^. This leads to the following truncation parameters 
Zj = ^ , and the truncated characters needed for the counting become 



r„(x;g,fc)=^g5(-( 



I -a) 



fc-1 

n 

i=l 



+ ((I,_i-C,_i)-a), 



(41) 



To obtain the symbols { ]? }fc which are needed for the counting, one has to take the limit q ^ 1 
of the prefactor of x^ in eq. (pl|). This results in 



fe-i 

E 

^ iai=F i=l 



n 



((ifc- 



1 — 'Ck-i] 
ai 



(42) 
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With this result, we arrive at the following counting formula for the Read-Rezayi states (for general 
k) 



with n = kAN^. To make the above (in particular the symbols {}k) more explicit, we will discuss 
the k = 2 (i.e. the pfaffian) and fc = 3 cases. For the pfafhan state counting, we need to know the 

n 

symbols {}2- Eq. ( p^ ) with k ~ 2 gives {f}2 — (|,)- Of course, this is just the result already 
found in . Note that our notation is slightly different with respect to the one used in ^ . In 
our notation, we denote the number of created quasiholes by n. In |l2| n denoted the number 
of extra fluxes, which is denoted by AA^^ in this paper. 

Although the method described above seems to be unnecessarily complicated to reproduce this 
result, it is very useful for obtaining closed expressions for k > 2. As an illustration, we will discuss 
the case k — 3, and compare our results with For fc = 3, the polynomials are given by the 
following expression 



F„(x;q,3) = 



X 

a^b 



rn a+2bl 




V2n 2a+b-\ 


3 3 




3 3 


a 




h 



(44) 



Indeed, these polynomials reduce to the ones in upon setting q = \. The symbols {ja are now 
easily written down 



n 



/ n _ a+2b \ / 2n _ 2a+b 
a+2b=F ^ ^ ^ 



(45) 



Note that only a finite number of terms contribute to the sum in eq. ([451). In fact, this is true for 
all the symbols ( ^2|) with n finite. 

The fusion rules for the spin field a which is part of the quasihole operator at level fc = 3 (see 
Q), can be encoded in a Bratteli diagram with the same structure as the diagram^ (note that the 
fields differ, of course). This is a consequence of the rank-level duality su(2)3 ^ su(3)2 (see §16.6 
in [0). Thus the total intrinsic degeneracy for the fc = 3 Read-Rezayi state with n quasiholes is 
given by dn = J-{n — 2). Indeed, by summing the symbols {f}3 over F, this result is reproduced. 

The angular momentum multiplets can be found in the same way as described in section |^. Let 
us note that for fc = 1 the only degeneracy factor remaining is (^^") , which is precisely the orbital 
factor for the Laughlin states with quasiholes present. This was of course to be expected, as the 
fc = 1 Read-Rezayi states are in fact the Laughlin states. 

To conclude the discussion on the counting for the Read-Rezayi states, we would like to mention 
that the numerical studies as presented for fc = 3 in |^ are in complete agreement with the counting 
formulas. At this point, no numerical results are available for fc > 4. In the following section we 
will turn our attention to the counting of the NASS states for general level. 

9. Counting formulas for the NASS states 

In this section, we describe the counting for the NASS states at general level fc. We will closely 
follow the procedure of the previous sections, that is, we start by writing down the chiral character 
corresponding to the 5u(3)fc/u(l)^ parafermions ^ 

E„i(a-Ck-i-a+b-Ck-i-b-a-Ck-i-b) 
n t \ f \ ^'i'''^2 ' . (46) 
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where we used the same notation as in eq. (|37|). This character is of the UCPF form with the 
K-matrix equal to the K-matrix of the su(3)fc/u(l)^ parafermions: K = ( V ) ® ^k~i [0- 
recursion relations corresponding to the basis of this theory can be written in the following way 



^(l,m) {xi,X2; q,k) 
Y{i,in){xi,X2;q,k) 



Y, 



i(k- 



(1— ei,m) 



(47) 



Once again, we solve the recursion relations by matching these recursions to eq. (^0|). The 
truncated characters take the form 



Y(i^m)ixi,X2;q,k) 



: ^ ^^(a-Ck-i-a+b-Ck-i-b-a-Ck-i-b) ^ 

ai.hj 
k-1 



n 

i=l 
k-1 

n 



-Cfe_i) -a+iCfe-i -b) 



(m+ (Ife_i 



Ik-i) ■ b + |Cfc_i ■ a)^. 



(48) 



We continue by specifying the parameters li and rrij. We have to use the same construction as 
in the RR case, with the difference that we now need the number of spin-up and down particles 
(denoted by rif and ni) created by the excess flux. Using k = ^ and mj = results in 



Y(^„j,„^)ixi,X2;q,k) 



q|(a-Ck-i-a+b-Ck-i-b-a-Ck-i-b) ^ 

ai,b-i 



k-1 



1=1 

k-1 



((Ife_i-Cfe_i)-a+iCfe_i-b). 



n 



fe 



(49) 



From eq. 

•^1 -^2 



((Ifc_i-Cfc_i)-b+iCfc_i-a)^ 

we obtain the symbols { 'p^^k by taking the limit 5 ^ 1 of the coefficient of 

^ + ((Ife_i-Cfe_i)-a+iCfe_i-b),^^ 



Fi F2 



k-1 

E n 

Y: iai=Fi i=l 
k-1 /jn-f 



n 

i=i 



+ ((I;,_l-C;,_l)-b + 

bj 



2^k-i ■ aj^. 



We know have specified all the ingredients of the counting formula for the NASS states 

#NASs(iV,A7V^,fc) - J2 



A''T,i:"T,i.-P"l.2 



Fi F2 



k ' "T 



N1-F2 

k 



(50) 



(51) 



where the prime on the sum indicates the constraints + = N, n-^ + rii ~ 2kAN^ and 
N-^ — Ni = III — The last constraint is a necessary condition for the state to be a spin-singlet 
(for more information on the constraints, see |To[). 

The case k = 2 was explicitly discussed in section 0. For fc = 1 only the orbital degeneracy 
factors remain, and we obtain the counting formula for a particular class of Halperin states p^ . 
Indeed, for k — 1, the NASS states reduce to the spin-singlet Halperin states. As already mentioned 
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in section |j, the counting formula (Bl|) with the symbols (^0|) exactly reproduces the results of the 
diagonalization studies for k — 2 [[10[. For fc > 3, no numerical results are available at the moment. 

10. Discussion 

In this paper, we explained in which way parafermions can be used to define clustered quantum 
Hall states. The statistics of these parafermions is needed to understand the energy spectrum of 
the clustered states in the presence of quasiholes, as obtained from a numerical diagonalization 
study. We obtained explicit formulas for the symbols {}k for two classes of clustered states, needed 
in the counting formulas. From these formulas, also the spin and angular momentum multiplet 
structure of the quasihole degeneracies can be extracted. 

One can say that the parafermion statistics is a crucial part in the understanding of the ground 
state properties of a system of interacting electrons (via an ultra local interaction) on the sphere in 
the presence of a magnetic field. In fact, we know of no other way in which this energy spectrum can 
be understood. The observation that ground state properties of a system of interacting electrons 
needs the knowledge of parafermion statistics is interesting by itself. 

Recently, another class of paired (clustered) spin-singlet states was proposed |^ . An interesting 
property of these states is that the fundamental excitations over these states show a separation 
of the spin and charge degrees of freedom. We believe it should be possible to repeat the present 
analysis for these newly proposed states, though it will be more difhcult, because the underlying 
parafermions are related to so (5), a non simply-laced Lie algebra. We leave this as a challenge for 
future work. 
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